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We report, for the first time, tlie generation of high-order harmonics in a spectral range between 200 eV and 
1 keV with an unusual spectral property: only every 4*^ (4i -h 1, i e K) harmonic line appears, whereas the 
usual high-harmonic spectra consist of every odd (2i -H 1) harmonic. We attribute this unique property to 
the quantum path interference of two extended electron trajectories that experience multiple re-scattering. 
In the well-established theory, electrons emitted via tunnel ionisation are accelerated by a laser field, return 
to the ion and recombine. The acceleration typically lasts for less than one optical cycle, and the electrons 
radiate in the extreme ultraviolet range at recombination. In contrast, for extended trajectories, electrons are 
accelerated over two or more optical cycles. Here, we demonstrate that two sets of trajectories dominate and 
provide substantial contributions to the generated soft X-ray radiation because they fulfil the resonance 
condition for X-ray parametric amplification. 

In the past several years, high-harmonic generation (HHG) has become a promising method of generating 
coherent radiation in the extreme ultraviolet (XUV) and soft X-ray (SXR) spectral ranges. HHG converts the 
frequency of intense laser pulses from the visible or near-infrared region into its harmonics whilst maintaining 
the spatial and temporal coherence of the driving laser pulses. The generated short-wavelength pulses have 
durations that range from femtoseconds down to a few tens of attoseconds^ and the spectra can span the 
wavelength range from 100 nm to far below 1 nm^"^. This type of radiation paves the way for high- resolution 
time- resolved spectroscopy in the keV spectral range^'^'^. 

Shortly after the experimental discovery of HHG^, a theoretical model was developed^ that is able to describe 
the experimental findings very accurately. According to this model, a strong but non-relativistic laser field 
generates the harmonic field from single atoms or molecules in the gas phase. The observed harmonic field is 
the coherent superposition of independently generated fields, and the macroscopic build-up is affected by the 
phase-matching conditions^'^°. The model described in Ref. 8 separates the HHG process into three distinct steps. 
First, in a strong laser field, an electron is released from an atom via tunnel ionisation. Next, the electron is 
accelerated in the laser field and gains energy within a time frame of approximately y4 optical cycles. Finally, the 
electron returns to the atom and recombines, and the excess kinetic energy is radiated as an XUV photon. This 
model relies on the strong-field approximation, i.e., the Coulomb field of the ion is neglected during the electron 
acceleration. Furthermore, because of the non-relativistic intensity, the magnetic field of the laser beam is also 
neglected. If the magnetic field is not neglected, one must consider a small, but non-negligible (even for non- 
relativistic intensities), acceleration of the electron in the direction of the laser propagation. The resulting 
longitudinal displacement of the re-colliding electron compared to the ion is referred to as the impact parameter 
b. 

The generated high-harmonic spectrum contains only the odd harmonics (q = 2i + 1, i e K) of the generating 
laser frequency with a well-defined cut-off (highest) energy. The preferred method of extending the cut-off into 
the soft X-ray region is to increase the intensity (I) or the wavelength (k) of the driving laser pulses, as the cut-off is 
proportional to X^I. However, for longer driving wavelengths, the electron wave packet spreads out farther, 
resulting in a rapid decrease in the recombination probability and a reduction in the high-harmonic generation 
efficiency^ On the other hand, an increase in the intensity or wavelength of the driving laser increases the 
impact parameter b (^ oc ny^. Again, the higher the impact parameter, the lower is the recombination 
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probability and, hence, the HHG conversion efficiency. However, it 
has been theoretically proposed^^ that for higher laser intensities, up 
to the relativistic regime, extended electron trajectories can overcome 
the limitation imposed by the impact parameter, and hence, harmon- 
ics can be generated up to the hard X-ray range; furthermore, the 
interference of the multiple re- scattered electrons originating from 
different optical cycles determines the temporal build-up of HHG^^. 
An electron moving along an extended trajectory returns to the ion 
much more than Va optical cycles after emission, which is the latest 
possible return time for conventional HHG. 

In this study, we present, for the first time, experimental evidence 
of HHG radiation from extended electron trajectories. The measured 
spectra, which extend to 1 keV, can be explained only by the pres- 
ence of X-ray parametric amplification^^'^^ (XPA) combined with 
HHG from electrons returning to the ion two or more optical cycles 
after ionisation. In our recent experiments, high-order harmonic 
radiation was generated using 30-fs-long laser pulses with peak 
intensities of up to 1 X 10^^ W/cm^ and a centre wavelength of 
800 nm. For this parameter range, many possible electron traject- 
ories exist. However, the resonance condition for XPA is best 
matched by two distinct sets of electron trajectories that produce 
two X-ray pulses in every optical half-cycle. The quantum path inter- 
ference between the two extended electron trajectories results in 
harmonic spectra that contain only every 4* harmonic (q = 4i + 
1) in the spectral range of a few hundred eV to keV. 

Results 

In this section, we present experimental results that demonstrate the 
unusual behaviour of the generated high- order harmonic spectra, 
and we then propose an explanation for this unexpected behaviour 
by invoking the combined action of two processes: (i) the quantum 
paths of two electron trajectories interfere during the generation 
process, and (ii) the two trajectories fulfil the conditions for X-ray 
parametric amplification. 

Measured uncommon HHG spectra. A series of spectral measure- 
ments was performed using the JETI Ti:sapphire laser system of the 
university in Jena, Germany. It delivers 30-fs-long pulses with 
energies of up to 1 J at a central wavelength of 800 nm and a 
repetition rate of 10 Hz. The high laser-pulse energy permitted a 
loose focusing geometry (f/100), and a peak intensity of up to 1 X 
10^^ W/cm^ in a heUum gas target was achieved. The HHG spectra 
were measured in the spectral range of 0.2-1 keV using an X-ray 
spectrograph (details can be found in the Methods section). A 
recorded X-ray spectrum is presented in Fig. 1(a) to permit the 
detailed examination of the spectrum in the range between the 
200* and 360* harmonic, corresponding to an energy range of 
approximately 300 to 600 eV, respectively. In this range, we 
observed primarily harmonic lines corresponding to the harmonic 
orders of q = 4i + 1 (where i is an integer). The harmonic lines are 
very pronounced in the range of the 20 to 22 orders and between 
the 269* and 289*. In the range between harmonic orders 221 and 
261 and near the high-energy cut-off, the signal is much weaker 
because of the lack of phase matching; this phenomenon will be 
discussed later. Above harmonic order 321, the visibility of the 
harmonic lines is weaker as a consequence of the decreasing 
spectral resolution of the spectrograph and the limited sampling 
rate. The small anomaly near harmonic order 293 is attributed to 
the L-edge absorption of the thin Ti-filter in the beam path. 

The full spectra of the generated high-order harmonic lines were 
further examined by calculating their Fourier transforms. The results 
of the transformation clearly reveal a quarter-optical-cycle peri- 
odicity for signals in the spectral range of up to 1 keV (1.2 nm) for 
a gas backing pressure of 0.5 bar, as shown in Fig. 1(c). The small 
deviation from the expected exact % -cycle periodicity originates 
from the blue shift of the harmonic lines, which is a consequence 



of the blue shift of the driving laser wavelength caused by the fast 
ionisation in the gas jet^^ and the atomic phase contribution^^. The 
additional frequency components in the wavelength range below 
3 nm (—267* harmonic) can be well explained as numerical artefacts 
caused by the low signal-to-noise ratio of the measurement in this 
spectral range. Furthermore, additional spectral modulation arises 
from the Ti foil in the beam path. The spectral range of interest 
extends well above the Ti L-absorption edge, so we expect a strong 
spectral modulation in the near- edge region (XANES) and the 
extended range above the edge (EXAFS)^°. In the results summarised 
in Fig. 1(c), we can also identify weaker modulations that correspond 
to periodicities of V2 and optical cycles (green dashed lines), which 
can be explained as higher- order Fourier components attributable to 
the non-sinusoidal shape of the harmonic lines, especially at longer 
wavelengths. Further components originate from conventional har- 
monics with a half- optical- cycle periodicity (blue dotted line) gen- 
erated in the low-intensity leading pedestal of the laser pulse. This 
conclusion is supported by the measurement presented in Fig. 1(b). 
Here, the peak intensity was reduced by more than one order of 
magnitude to approximately 5 X 10^^ W/cm^, and we observed every 
odd harmonic in the low-photon-energy range, as expected for con- 
ventional HHG. The additional lines that deviate from the predicted 
positions and the line at the position of the 249* harmonic in 
Fig. 1(a) are attributed to atomic/plasma emission lines. 

X-ray parametric amplification. Here, we briefly summarise the 
important features of the X-ray parametric amplification (XPA) 
process that contribute to the present experimental findings. 

In our model (described in detail in Ref 16 and 17) based on the 
conventional three-step model of high-order harmonic generation^'^, 
the Coulomb field of the ion plays an important role in the third step, 
i.e., at the instant when the electron re-collides (re-scatters) and the 
X-ray photon is generated. In the equation that describes the inter- 
action between the electron and the ion, the Coulomb force is 
expanded into a Taylor series, and the contributions of the different 
orders can be discussed separately. The 0*-order term defines a 
resonance frequency co^ ~e^/(47r£o/tZe^^), where b is the impact para- 
meter, i.e., the distance of the scattering electron from the ion. The 
first-order term describes a parametric interaction, i.e., a high-order 
difference-frequency generation, namely, oj^ = 2qwi — Wq, where Wi 
and cOq are the angular frequencies of the first-order Ei (the driving 
laser) and the q*- order harmonic Eq field. This parametric coupling 
between the laser and X-ray fields appears in the form of 
dzEqCcEf^E*, giving a contribution to the signal only if an X-ray field 
Eq is already present, e.g., from conventional HHG. The impact 
parameter also plays a very important role in the X-ray parametric 
amplification process; if it is of the correct order of magnitude, the 
conditions for parametric resonance are fulfilled, resulting in a sub- 
stantial amplification of the short-wavelength signal generated via 
conventional HHG. This amplification process has been experiment- 
ally observed and theoretically described in Ref 16 and 17, respect- 
ively, and is known as X-ray parametric amplification (XPA). 

Extended electron trajectories. In this section, we describe how 
extended electron trajectories can be created, how they behave and 
how they can support XPA in the few-lOO-eV range by examining 
the distance of the released and accelerated electron from the ion. 

In a strong laser field, atoms are ionised via tunnel or above-barrier 
ionisation. Subsequently, the freed electrons are accelerated in the 
laser field. The probability of tunnelling ionisation depends on the 
electric field, so most electrons are released near the maximum of 
the electric field. The motion of the freed electrons in the laser field 
has been calculated by considering the electric and magnetic fields of 
the laser beam and the Coulomb potential of the ion^^'^^ (for details, 
see the Methods section). Note that in the context of XPA, it is not 
necessary to consider a softened Coulomb potential to correctly 
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Figure 1 | Measured spectra with unusual spectral lines, (a) Only harmonic lines of order q = 4i + 1 appear in the measured HHG spectra in the spectral 
range of 300-600 eV for a laser peak intensity of 8 X 10^^ W/cm^. (b) At a peak intensity of approximately 5X10^"^ W/cm^, the harmonic lines appear at 
every odd harmonic position, (c) The harmonic lines and their quarter- optical- cycle periodicity can be resolved up to ~1 keV (1.2 nm) at 0.5 bar 
and 8 X 10^^ W/cm^. High-order Fourier components (highlighted with green dashed lines) also appear. Standard harmonics of half- optical- cycle 
periodicity (highlighted by blue dashed lines) are generated in regions with a lower intensity in space and/or time. 



predict the experimental results. In our model, the electron traject- 
ories are calculated classically, as in many high-harmonic generation 
theories^^. At re-collision, when the electron is described quantum 
mechanically, the impact parameter represents the distance to the 
centre of gravity of the electron wave packet and can be easily 
included in the calculation for the atomic polarisation^^. 

To match our experimental conditions, we considered a 30-fs-long 
laser pulse with a central wavelength of 800 nm and a peak intensity 
of 5 X 10^^ W/cm^orS X 10^^ W/cm^ Furthermore, the calculations 
were performed at the leading edge of the laser pulse, where atoms 
are usually ionised and harmonics are generated. In Fig. 2 (a) -2(d), 
we show the distance of the electron from the ion during acceleration 



in the electric and magnetic fields of the laser pulse. The details of the 
model used and the manner in which the initial distance and speed of 
the electron were determined are described in the Methods section. 
For a laser intensity of 5 X 10^^ W/cm^ (Fig. 2(a) and 2(b)), we show 
the trajectories for electrons emitted 6 or 7 optical half cycles (OHC) 
before the intensity maximum, whereas for a peak intensity of 8 X 
10^5 W/cm^, we considered electron trajectories originating 10 and 
1 1 OHC before the intensity maximum. These trajectories have been 
confirmed to support XPA. Furthermore, these electron trajectories 
were followed throughout several optical cycles, including several re- 
scattering events on the ion, where the re-scatterings were considered 
to be elastic. These electron trajectories must begin within a 20-as 
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Figure 2 | Extended electron trajectories. For He gas interacting with a 30-fs-long laser pulse with a central wavelength of 800 nm, extended electron 
trajectories were calculated for two intensities: (a), (b) 5 X 10^^ W/cm^ and (c), (d) 8 X 10^^ W/cm^. (a)-(d) The distance of the electron from the ion was 
calculated throughout several re-scattering events and plotted as a function of time, assuming that ionisation occurred a few optical half-cycles (OHC) 
before the intensity maximum. The scale of the colour bars is proportional to the probability of finding the electron at a certain distance from the atom. 
(e)-(h) At every re-scattering, the moment of re-scattering (re-collision) and the kinetic energy of the electron were also calculated and exhibit strong 
correlation. In every OHC, two re-collisions occur, and two soft X-ray pulses are generated. The scale of the colour bars is proportional to the 
probability of the electron having a certain kinetic energy at the time of collision, (i), (j), (1), (m) For both intensities, two corresponding electron 
trajectories with a re-collision time of approximately 2 optical cycle are presented; they fulfil the condition of X-ray parametric amplification in the 
corresponding X-ray energy range because they collide with suitably small impact parameters, (k), (n) The temporal evolution of the distance of the 
electron from the ion is also plotted separately for both chosen trajectories. 



time window (beginning 10 as, 15 as, 10 as and 5 as after the electric- 
field maxima of 6, 7, 10 and 11 OHC, respectively); otherwise, the 
electrons will not undergo multiple re-scattering, or, even worse, may 
not re-scatter at all, implying that neither HHG or XPA can take 
place. 

For all trajectories originating in the stated time window, the 
impact parameters at the first (bi) and second (ba) returns are too 
large for a high recombination probability. However, the electrons 
will interact with the ion via elastic scattering. So, with a high prob- 
ability, the electrons are further accelerated in the laser field, travel 
along extended trajectories, and return to the ion several more times. 
The impact parameter after the third (bs) or fourth (b4) return can be 
suitably small, not only enabling recombination but also fulfilling the 
condition for X-ray parametric resonance in the few-lOO-eV range. 
The resonance frequency Wr for XPA is a function of the impact 
parameter and can be calculated^^'^^ as CL)J~e^/(47r£oWe^^), where e 
and me are the charge and mass of the electron, respectively. The 
impact-parameter range for resonance is highlighted in yellow in 
Fig. 2(a) -2(d). When the electron trajectories pass through this 
range, then the condition for XPA resonance is fulfilled. This can 
occur two, three or possibly more optical cycles after ionisation. 
Furthermore, the electron must have gained sufficient kinetic energy 
at the instant of re-collision. The acquired electron energies are plot- 
ted in Fig. 2(e)-2(h), where the time windows for re-collision after 2 
and 3 optical cycles after ionisation are indicated. The re -collision 
energies have characteristic V shapes, implying one collision in every 
optical half cycle at low X-ray energies and two collisions at higher X- 
ray energies: two collisions in the second half of the ranges (Fig. 2(e) 



and 2(g), highlighted by grey and pink backgrounds) and another 
two collisions in the first half of the ranges (Fig. 2(f) and 2(h), high- 
lighted by grey and pink backgrounds). In summary, there are up to 
four collisions in every optical cycle. Within a certain energy range, 
between 200 and 400 eV in Fig. 2(e) and 2(f) and between 300 and 
500 eV in Fig. 2(g) and 2(h), these collisions are essentially uniformly 
distributed, leading to the full destructive interference of every sec- 
ond harmonic line, thereby yielding spectra in which only every 4*^ 
harmonic line is present and explaining the measurement results. 

The shapes of the two corresponding electron trajectories, with re- 
collision times of approximately 15/8 and 17/8 optical cycles, are 
shown in more detail in Fig. 2(i) and 2(j), respectively, for a laser 
intensity of 5 X 10^^ W/cm^ and in Fig. 2(1) and 2(m), respectively, 
for a laser intensity of 8 X 10^^ W/cm^. The corresponding distances 
of the electron from the ion are summarised in Fig. 2(k) and 2(n) for 
the two different intensities. For both laser intensities, the corres- 
ponding trajectories appear very similar; only the instants of ionisa- 
tion (18 as and 7 as or 33 as and 13 as, respectively, after the local 
electric-field maximum) are somewhat different. As a consequence 
of the quantum path interference between these two electron traject- 
ories, the generated harmonic spectrum will consist only of harmon- 
ics of the order q = 4i + 1, where q is the harmonic number, in 
contrast to the conventional high-harmonic spectrum. 

At lower laser intensities, e.g., 5 X 10^^ W/cm^, the kinetic energies 
of the electrons are sufficient for generating HHG radiation in the 
spectral range below 150 eV. Furthermore, at such a moderate laser 
intensity, the impact parameters at the first returns (bj are very 
small. Therefore, if the electrons do not recombine, then their tra- 
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jectories will be strongly perturbed by the Coulomb potential of the 
ion, circumventing further returns; consequently, extended electron 
trajectories with multiple re-scattering do not exist. Nevertheless, in 
this low-intensity regime, high-order harmonics are generated for 
electrons that return —7/8 optical cycles after ionisation, which are 
commonly known long trajectories, and the corresponding impact 
parameters are in a range that fulfils the condition for XPA in the 
spectral range below 150 eV. Here, ionisation and recombination 
can occur only once in every optical half-cycle, and the observed 
high-harmonic spectrum consists of odd harmonics (q = 2i + 1), 
in full agreement with the conventional high-harmonic spectrum 
(see Fig. 1(b)). 

Phase matching and XPA. In Fig. 1(c), the low-frequency Fourier 
components were filtered out so that they would not overwhelm the 
important high-frequency components. However, the lower 
components carry important information regarding the phase- 
matching condition in HHG. The phase mismatch modulates the 
amplitude of the harmonics following a sinc^((p(p)/2) function, 
where (p(p) = 2npQ^qp; q is the harmonic order, p is the gas 
pressure and p^^ is a phase-matching parameter. Corresponding 
sinc^ envelopes are plotted in Fig. 1(a) (dotted black line) along 
with the transmission of the thin-film Ti filter used in this study, 
and the sin^ envelopes are also plotted in Fig. 3(a), which presents 
several spectra from the series recorded at various pressures. The 
Fourier transforms of the spectra contain two distinct peaks at very 
small fractions of an optical cycle (see three examples in Fig. 3(b)). 
The peak positions from every evaluated spectrum are plotted in 
Fig. 3(c), and this figure exhibits a linear dependence for both 
spectral contributions, as expected from theory. The blue line 
crosses the axis near zero, as can be expected for XPA (see Eq. (32) 
in Ref 17). The red peaks should originate from the conventional 
HHG contribution in the generated spectra. The slope of the fitted 
line for the XPA contribution is smaller, which indicates a better 
phase-matching condition; and a p^^ phase-matching parameter 
of po = 64 ± 6 bar was determined from this slope. Using this 
parameter, the calculated pressure dependence of the XUV signal 
in the presence of XPA was fitted to the measured data (see Fig. 4) in 
various spectral ranges. For the fitting, the electric field of the q* 
harmonic as a function of the pressure was estimated to be^^ 

CP 



Eq(p) = aqexp(Gq(p)) 



exp (i(p(p') — Gq(p'))dp\ with the pressure- 



cP' 

dependent gain Gq(p') exp(—p^/p^^)dp, and aq denotes the 

2 Jo 

signal amplitude of the phase-matched high harmonics without 
XPA. The measured and fitted curves in three different photon 
energy ranges exhibit perfect agreement between the measurement 
and the prediction of our XPA model, as shown in Fig. 4. After 
estimating a saturation pressure of ps = 0.7 bar from our 
measurement, only two free fitting parameters remain in our 
model, namely, the phase-matched HHG intensity (oca^) and the 
XPA gain (gqo). The parameters obtained for the best fits in various 
photon -energy ranges are summarised in the inset of Fig. 4. The 
measured gain parameter is nearly independent of the photon 
energy; this result is consistent with our model, assuming a 
uniform electron distribution at the position of re- collision 



Discussion 

From the presented measurements, we can draw the following major 
conclusions. For laser intensities in the range of 10^^-10^^ W/cm^, 
high-order harmonics can be generated, and the X-ray yield is 
enhanced via parametric amplification from extended electron tra- 
jectories. For these extended trajectories, the time between ionisation 
and recombination is approximately two or more optical cycles, and 
in every optical cycle, four X-ray pulses are generated. The prolonged 
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Figure 3 | Phase matching, (a) In the measured spectra, long-period 
intensity modulations can be observed (black dashed lines) with pressure- 
dependent periods. This modulation is identified as a consequence of the 
non-phase-matched generation of the high harmonics, (b) The long- 
period modulations produce two distinct peaks in the Fourier transforms 
of the spectra at positions that depend on the gas pressure, (c) In 
accordance with the pressure dependence of the peak positions, an XPA 
contribution and a non-phase matched HHG contribution can be 
recognised in the generated spectra. 



time spent as a free electron results in a more pronounced spreading 
of the electron wave packet^^ thus lowering the recombination prob- 
ability. However, the spreading of electron wave packet can be com- 
pensated by Coulomb focusing^\ 

The findings presented in this study provide an improved under- 
standing of high-order harmonic generation and X-ray parametric 
amplification at higher laser intensities and will allow the establish- 
ment of guidelines for extending the XUV signal to higher photon 
energies with increased efficiency. It is therefore worth mentioning 
that our X-ray source, at approximately 1 keV, provides a photon 
flux comparable to those of other recently published approaches 
using driving lasers with longer wavelengths^. Such an efficient 
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source paves the way for attosecond science far beyond 100-eV 
photon energies and/or using shorter laser wavelengths for improved 
efficiency and photon yield^^'^^. 

Methods 

Electron-trajectory calculation. In a strong laser field, atoms are ionised, and the 
released electrons are accelerated in the laser field. The motion of the electron can be 
described by the classical equations of motion in the static Coulomb field of the ion 
and the Ei electric and B i magnetic fields of the driving laser as follows: 



where r is the distance of the electron from the ion. For the calculation, we consider a 
linearly polarised laser beam propagating in the z direction and polarised in the x 
direction. Thus, the electric field has only an x component, Ei = (£i, 0, 0), and the 
magnetic field has only a y component, Bi = (0, 5i, 0); therefore, the motion of the 
electron remains in the x-z plane. The trajectory is defined by the initial conditions 
{ti,Xi,Vi}, where ti is the moment of ionisation and Xj and Vi are the position and 
velocity, respectively, of the electron at tj. In the calculation, we consider laser pulses 
with a central wavelength of = 800 nm, a Gaussian temporal shape with a duration 
of AT = 30 fs (FWHM) and a peak intensity of 5 X 10^' W/cm' or 8 X lO^^ W/cm\ 
in accordance with the experimental conditions. Driven by the strong electric field, 
the atom can be either tunnel ionised at lower intensities (see Fig. (5a)) or can be 
ionised above the barrier (see Fig. 5(b)). Because of the field dependence of the 
ionisation probability^ most electrons are released near the field maximum in every 
optical half-cycle (OHC). At such high intensities, the ionisation rate is very high, so 
almost all atoms are ionised by the leading edge of the laser pulse within a few optical 
cycles prior to the intensity maximum. The time of ionisation is defined as 



= OHC-— +fo, 
2c 



(2) 



where to is measured from the maximum/minimum of the electric field, and at this 
time, the electric field is 



Ei = Eq exp 



-2ln(2)(A) 



2kc 
cos I — ti 



(3) 



For our estimation, we consider helium atoms. The energy required to remove the 
first electron from the He atom is 24.6 eV, and the energy required to ionise the He"^ 
ion is 51.4 eV, i.e., therefore 76 eV is required to fully ionise a He atom. This means 
that the two electrons in the He atom in the 1 s^ ground state are at the — 38 eV level 
(see Fig. 5(a) and 5(b)). For tunnel ionisation (Fig. 5(a)), one of the electron moves 
down to the —51.4 eV ground state of the He"^ ion, whilst the other moves up to the 
— 24.6 eV binding energy ( — Wb) and is released via tunnelling through the potential 
barrier at the level of the binding energy. This can be written as 



47reo|xf| 



- eEiXi ■■ 



(4) 



and provides a solution for the location Xi, outside the barrier, where the electron 
appears after tunnelling: 



Wb/e 
2Ei 



e\Ei\ 



nso(Wb/ey 



(5) 



According to the principle of energy conservation, the electron released at Xi must have 
zero velocity, Vi = 0. The electron can tunnel at a lower energy level (thin arrows in 
Fig. 5(a)), and the remaining energy will be the kinetic energy of the electron with non- 
zero initial velocity Vj' > 0. However, the probability of ionisation through a broader 
barrier is much smaller, so ionisation at the binding-energy level should dominate. 
When the electric field satisfies the condition 



|£il> 



(6) 



4n£Q 



— r-e(Ei+rxBi), 



which is 104.9 V/nm for He, the atoms are ionised over the barrier, and Xj is defined at 
the maximum of the barrier (see Fig. 5(b)): 




Figure 5 | Tunnel and above-barrier ionisation of He atoms, (a) For an electric-field strength of lower than 104.9 V/nm, electrons (blue circles) from the 
ground state can leave the He atom via tunnel ionisation through the potential barrier (pink). The remaining electron of the He^ ion will then be 
much more tightly bound, (b) For higher electric-field strengths, the electrons will be freed above the barrier. After ionisation, the electrons will initially 
have the position Xi and the velocity Vj. 
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-sign(Ei) 



4:K£o\Ei\' 



(7) 



Furthermore, the remaining energy is released as kinetic energy and defines the initial 
velocity: 



-sign(Ei)^ 



2Wt 



(8) 



For our experimental conditions, the atoms were primarily ionised via above-barrier 
ionisation, so Eq. (7) and (8) define the initial distance and velocity. After ionisation, 
the electron is accelerated by the laser field and passes by the atom at a certain distance 
one or two times during every optical cycle. This distance is called the impact para- 
meter b, which is relatively small, so the bypass is referred to as re-collision. In the 
calculations, we followed the trajectory of the electron throughout several optical 
cycles and re-collisions. A set of trajectories was calculated for various ionisation 
times, with a resolution of 20 trajectories per attosecond. In Fig. 2, the distance of the 
electron at every time of ionisation (to) and the kinetic energy of the electron at every 
re-collision are plotted. 

HHG source. The gas target consisted of a Ni tube with an interaction length of 1 mm 
backed with helium in a pressure range of 0. 1 to 1 .2 bar. To reduce the gas load and the 
reabsorption, the target was operated with a pulsed valve at 10 Hz. The experimental 
conditions — He as the target material, the high ionisation rate that corresponds to a 
laser intensity of 8 X 10^^ W/cm^ and a gas-jet position that ensured that the atomic 
phase and the Gouy phase would compensate each other — were carefully chosen and 
fully ruled out the possibility of pressure-induced phase matching. In certain spectral 
ranges, this phase-matching scheme exhibited a similar scaling as XPA (Ref 23), 
suggesting the possibility of a different explanation of the experimental observations. 

HHG detection. After the fundamental laser light was blocked with thin Al (200 nm) 
and Ti (100 nm) foils, the HHG spectra were measured using an X-ray scanning 
monochromator (McPherson 248/310 G) equipped with a 1200 -groove/mm grating 
and a photo multiplier (Channeltron 4751 G). The photomultiplier was operated at a 
supply voltage of —2.5 kV (SRS PS350) and read out by a lock- in amplifier (Signal 
Recovery 5209) with an integration time of 0.1 ms. Using this equipment, we were 
able to record harmonic spectra in the spectral range of 0.2-2 keV with a resolution 
sufficient for the detection of single harmonic lines below 1 keV. 

Number of detected photons. The efficiency or sensitivity is not very precisely 
known for all components involved in the detection of the X-ray photons. Therefore, 
we opted for a device-independent method of estimating the detected number of X- 
ray photons. For our evaluation, we need only assume that Poisson statistics correctly 
describe the incoming photons. Then, we can estimate the number of photons N from 
the mean value (mean) and the standard deviation (sdev) using the following 
equation: N = (mean/sdev)^. For the estimation of N, we selected four well-resolved 
harmonic lines between 300 and 350 eV (see the spectrum in Fig. 1(a)) for a backing 
pressure of 0.3 bar. The signals within the harmonic lines were integrated, and the 
mean and sdev were calculated. Taking into account that every measurement point 
was produced by two laser shots, we obtained a value of 43 photons per harmonic line 
per pulse, and this value was used to correct the scales of Fig. 1(a) and Fig. 4. In the 
saturated regime, for a backing pressure above 0.8 bar (see Fig. 4), the number of 
photons was approximately 50 times higher, i.e., we detected approximately 2000 
photons per harmonic line per pulse. It should be noted that this method tends to 
underestimate the true number of detected photons because it assumes a driving laser 
with no fluctuations. 

Number of generated photons. To estimate the number of X-ray photons generated 
in the 300-350 eV range, we corrected the number of detected photons using the 
specifications of the elements in the beam path. According to the manufacturer, the 
photomultiplier used in this study has a quantum efficiency of approximately 10% in 
a spectral range of up to 1 keV. In the literature, we found the diffraction efficiency of 
the grating, which is on the order of 10"^ to 10"^. The transmission of the thin-foil X- 
ray filters (200 nm Al + 100 nm Ti) is approximately 10%, assuming no additional 
contamination on the filters. The size of the X-ray beam at the entrance of the X-ray 
spectrograph is unknown, making it impossible to correctly determine the fraction of 
photons transmitted through the input slit; therefore, we assumed a ratio between 0.1 
and 1 . The reabsorption of the generated X-rays in the gas jet and gas background was 
neglected, as it is relatively small in this spectral rage. All these contributions amount 
to an attenuation factor of between 10"^ and 10^. Therefore, in the 300-350 eV range, 
approximately 10^-10^ photons per line per pulse were generated. At 1 keV, the signal 
was approximately 2 orders of magnitude smaller, but the transmission of the filters is 
also higher, so we estimate that 10^-10*^ X-ray photons per line per pulse were 
generated. 



1 . Brabec, T. & Krausz, F. Intense few-cycle laser fields: Frontiers of nonlinear optics. 
Rev. Mod. Phys. 72, 545-591 (2000). 

2. Seres, J. et al. Source of coherent kiloelectronvolt X-rays. Nature 433, 596-596 
(2005). 

3. Seres, E., Seres, J. &: Spielmann, C. X-ray absorption spectroscopy in the keV range 
with laser generated high harmonic radiation. Appl. Phys. Lett. 89, 181919 (2006). 

4. Popmintchev, T. et al. Bright Coherent Ultrahigh Harmonics in the keV X-ray 
Regime from Mid- Infrared Femtosecond Lasers. Science 336, 1287-1291 (2012). 

5. Seres, E. & Spielmann, C. Time-resolved optical pump X-ray absorption probe 
spectroscopy in the range up to 1 keV with 20 fs resolution. /. Mod. Opt. 55, 
2643-2651 (2008). 

6. Seres, E., Seres, J. & Spielmann, C. Time resolved spectroscopy with femtosecond 
soft-x-ray pulses. Appl. Phys. A 96, 43-50 (2009). 

7. McPherson, A. et al. Studies of multiphoton production of vacuum- ultraviolet 
radiation in the rare gases. /. Opt. Soc. Am. B 4, 595-601 (1987). 

8. Lewenstein, M., Balcou, P., Ivanov, M. Y., L'Huillier, A. & Corkum, P. B. Theory of 
high-harmonic generation by low-frequency laser fields. Phys. Rev. A 49, 
2117-2132 (1994). 

9. Lompre, L. A. et al. High-order harmonic generation in xenon: intensity and 
propagation effects. /. Opt. Soc. Am. B 7, 754-761 (1990). 

10. Popmintchev, T. etal. Phase matching of high harmonic generation in the soft and 
hard X-ray regions of the spectrum. Proc. Nat. Acad. Sci. USA 106, 10516-10521 
(2009). 

11. Shiner, A. D. et al. Wavelength Scaling of High Harmonic Generation Efficiency. 
Phys. Rev. Lett. 103, 073902 (2009). 

12. Schiessl, K., Ishikawa, K. L., Persson, E. & Burgdorfer, J. Quantum Path 
Interference in the Wavelength Dependence of High-Harmonic Generation. 
Phys. Rev. Lett. 99, 253903 (2007). 

13. Kulander, K. C, Cooper, J. & Schafer, K. J. Laser-assisted inelastic rescattering 
during above-threshold ionization. Phys. Rev. A 51, 561-568 (1995). 

14. Hatsagortsyan, K. Z., Klaiber, M., Miiller, C, Kohler, M. C. & Keitel, C. H. Laser- 
driven relativistic recoUisions. /. Opt. Soc. Am. B 25, B92 (2008). 

15. Hernandez- Garcia, C. et al. Zeptosecond High Harmonic keV X-Ray Waveforms 
Driven by Midinfrared Laser Pulses. Phys. Rev. Lett. Ill, 033002 (2013). 

16. Seres, J. et al. Laser-driven amplification of soft X-rays by parametric stimulated 
emission in neutral gases. Nature Phys. 6, 455-461 (2010). 

17. Seres, J., Seres, E. & Spielmann, C. Classical model of strong-field parametric 
amplification of soft x rays. Phys. Rev. A 86, 013822 (2012). 

18. Rae, S. C, Burnett, K. & Cooper, J. Generation and propagation of high-order 
harmonics in a rapidly ionizing medium. Phys. Rev. A 50, 3438-3446 (1994). 

19. Shin, H. J. et al. Nonadiabatic blueshift of high-order harmonics from Ar and Ne 
atoms in an intense femtosecond laser field. Phys. Rev. A 63, 053407 (2001). 

20. Bressler, C. & Chergui, M. Ultrafast X-ray Absorption Spectroscopy. Chem. Rev. 
104, 1781-1812 (2004). 

21. Brabec, T., Ivanov, M. Y. & Corkum, P. B. Coulomb focusing in intense field 
atomic processes. Phys. Rev. A 54, R2551-R2554 (1996). 

22. Falcao-Filho, E. L., Gkortsas, V. M., Gordon, A. & Kartner, F. X. Analytic scaling 
analysis of high harmonic generation conversion efficiency. Optics Express 17, 
11217-11229 (2009). 

23. Kazamias, S. et al. Pressure-induced phase matching in high-order harmonic 
generation. Phys. Rev. A 83, 063405 (2011). 



Acknowledgments 

This study has been sponsored by the DFG grants TR18 PIG, TR18 P12, and SE 1911/1-1; 
the TMBWK grant B 715-08008 and the European Fund for Regional Development (EFRE). 
E.S. acknowledges support from the FSU grant "ProChance 2009 Al". The authors 
acknowledge the contributions of the JETI laser team. 

Author contributions 

J.S., E.S., T.K. and C.S. designed the experiments and wrote the manuscript; J.S., E.S., B.L., 
B.E. and B.A. conducted the X-ray experiments; J.S. developed the theory. All authors 
analysed the data and contributed to the completion of the manuscript. 

Additional information 

Competing financial interests: The authors declare no competing financial interests. 

How to cite this article: Seres, J. et al. High-harmonic generation and parametric 
amplification in the soft X-rays from extended electron trajectories. Sci. Rep. 4, 4234; 
DOI:10.1038/srep04234 (2014). 



l(cc^(D®(§)l ^^^^ ^o^^ is licensed under a Creative Commons Attribution- 

i^~j^:i«!r.».-i.J NonCommercial-ShareMke 3.0 Unported license. To view a copy of this license, 

visit http://creativecommons.Org/licenses/by-nc-sa/3.0 



SCIENTIFIC REPORTS | 4 : 4234 | DOI: 1 0.1 038/srep04234 



7 



